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Nitrogen, as the seventh most abundant element in the universe, 
is an important constituent of the atmospheres and interiors of 
planets such as the Earth and the surfaces of moons such as 
Triton. The phase diagram and equation of state of dense ni- 
trogen is therefore of interest in understanding the fundamental 
physics and chemistry of planetary processes and in discovering 
new materials. We predict stable phases of nitrogen at multi- 
TPa pressures, including a PA/nbm structure consisting of par- 
tially charged N^^ pairs and tetrahedra, which is stable in 
the range 2.5-7.1 TPa. This is followed by a modulated Fdd2 
structure at 7.1-11.5 TPa, which also exhibits significant charge 
transfer. The P4/nbm metallic nitrogen salt and Fdd2 modulated 
structure exhibit strongly ionic features and charge density distor- 
tions, which is unexpected in an element at such high pressures 
and could represent a new class of nitrogen materials. 

New materials under extreme conditions | Density functional theory | ab 
initio I Crystal structure predictions | polymerization of molecular solids 

The polymerization of molecular nitrogen under pressure has been 
very actively researched over the past two decades and has stim- 
ulated many experimental |T]|2l34 56|71[8l|9l and theoretical 
fra[TT1ll2llT3][T4llT5]|T6l[T7l[T8hl9.,20..2Tll2ai23^ 

investigations. Density functional theory (DFT) studies suggested 
that dissociation of nitrogen could occur at high pressures which 
were, nevertheless, attainable in diamond anvil cell (DAC) experi- 
ments. In a ground-breaking paper, Mailhiot et al. |13| predicted 
polymerization of nitrogen molecules under pressure, leading to the 
formation of the "cubic gauche" framework structure. After consid- 
erable efforts UlIU, cubic gauche nitrogen was finally synthesized 
by Eremets et al. (5), a decade after its prediction. The polymeric 
phase is quench recoverable to ambient conditions, where it exists 
as a rather unstable high-energy-density material (HEDM) |5|. Un- 
der ambient conditions an energy of roughly 1.5 eV/atom (~ 10330 
J/g) is released on converting cubic gauche nitrogen to molecular N2, 
which is more than twice the energy density of TNT (about 4686 J/g). 
Polymeric nitrogen |5 1, or high-N content salts |28 1, are therefore ex- 
pected to be excellent candidates as high-energy-density materials. 

The N=N triple bond is one of the strongest known and break- 
ing it requires surmounting a substantial energetic barrier. However, 
at high temperatures the triple bond breaks at pressures above about 
110 GPa 1 5 1, a much lower pressure than predicted for the single bond 
in H2 1 29 1 and the double bond in O2, which is predicted to survive 
up to 2 TPa 1 30 3 1 1. Once the N=N triple bond is broken, a wide va- 
riety of structures can be adopted, similar to phosphorus and arsenic, 
and many studies aimed at elucidating the phase diagram of nitrogen 
have been published. 

The set of candidate structures includes simple cubic, the iso- 
electronic "black phosphorus" and a-arsenic or cis-trans and Cmcm 
chain structures, etc. [ll[l3l[II][Il[l5l[l6l[l6l[l7l[Ill|24l Never- 
theless, only the "black phosphorus" structure seemed feasible, as 
it was predicted to become more stable than cubic gauche above 
~210 GPa 1131 [161 l22l . The advent of structure searching using 
DFT methods marked a significant advance; under conditions in 
which standard chemical intuition may no longer be reliable, new sta- 
ble polymeric nitrogen phases beyond cubic gauche have been pro- 



posed, including the layered Pba2 and framework P2i2i2i struc- 
tures [26, 25J. "Black phosphorus" was subsequently ruled out on 
energetic grounds. These phases have in common a large DFT band 
gap that hardly closes with pressure. Nitrogen compounds could form 
metallic phases at high pressures and/or temperatures, which could 
have important consequences for understanding planetary processes. 

Large compressions normally lead to an increase in the atomic 
coordination numbers and metallization, as predicted in, for exam- 
ple, hydrogen |29| |32|, carbon |33 34|, and oxygen |30|, as well 
as observed in many other systems. |35| The coordination number 
of the atoms in solid nitrogen increases from one to three on trans- 
formation from the molecular to cubic gauche forms and three-fold 
coordination persists in higher pressure structures. 

All high-pressure nitrogen structures proposed so far are insulat- 
ing. Although an increase in coordination number with pressure is 
physically reasonable, it is by no means universal. For example, alu- 
minum transforms into more open structures at TPa pressures as the 
valence electrons move away from the ions in the formation of "elec- 
tride structures" |36| although, as must occur in a pressure-induced 
equilibrium phase transition, the density increases at the transition, 
and well-packed structures are strongly disfavored in oxygen up to 
at least 25 TPa |30|. The coordination number of nitrogen has not 
been determined at pressures above 400 GPa, and a metallic phase of 
nitrogen has not yet been reported. 

Materials under TPa pressures are important in planetary science, 
for example, the pressure at the center of Jupiter is estimated to be 
about 7 TPa |37|. Dynamical shock wave |37 38, 39| and ramped 
compression experiments I40II41I are increasingly being used to in- 
vestigate materials at TPa pressures. Even more extreme conditions 
are attainable in laser ignition experiments 1421 . Experimental de- 
terminations of structures are not, however, currently possible at TPa 
pressures and therefore theoretical predictions are particularly im- 
portant. I33ll32|[36ll34ll30l In this work we focus on the structural 
phases of nitrogen at multi-TPa pressures. 

Results 

Crystal structures and charge transfer. We have used ah initio ran- 
dom structure searching (AIRSS) (43' '441 and DFT methods to find 
candidate structures of nitrogen up to multi-TPa pressures. These 
searches have enabled us to identify novel candidate nitrogen struc- 
tures that are more stable than previously known ones. Four of them 
are thermodynamically stable within certain pressure ranges, namely. 
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PA/nbm (Z=14), Fdd2 {Z=4S), R3m (Z=l) and Mi/amd {Z=4). 
Most interestingly, the tetragonal Pi/nbm structure is characterized 
by the presence of N2 pairs and N5 tetrahedra, while the orthorhom- 
bic Fdd2 structure is modulated, as depicted in Fig.[TJa) and (b). 

At 2.5 TPa, the N-N bond length within the N5 tetrahedra of 
P4:/nbm is 1.13 A, and the shortest distance between the corners of 
the tetrahedra and N2 dimers is 1.26 A. The N-N bond length within 
the dimer of about 1.17 A is significantly shorter than the N-N sep- 
aration between adjacent dimers (about 1.29 A), which implies that 
the dimers are well separated. 

The unique structure of the PA/nbm phase led us to investi- 
gate the charges on the different atomic sites. We used the wave- 
function-based MuUiken population and charge-density based Bader 
1 45 1 methods to analyze the charge transfer. As summarized in Table 
[T] the two methods lead to the same conclusions; the N5 tetrahedra 
are negatively charged and serve as anions while the N2 dimers are 
positively charged and act as cations, and thus the PA/nbm structure 
resembles an all-nitrogen salt. The resultant charges (about ±0.22e 
from Mulliken analysis and about ±0.37e from Bader analysis) are 
substantial, although they depend somewhat on the definition of the 
atomic charges. 

The complicated Fdd2 structure has a 48-atoms conventional 
cell, and the atoms are arranged in undulating layers, as depicted in 
Fig. [TJb). The charge distribution of Fdd2 also exhibits large dis- 
tortions. As shown in Fig.[T|b) and Table [T] although all the atoms 
are at 16b sites, each orbit has a different character: there is a group 
of strong electron acceptors (Bader charge (5s=-0.36e, red), a group 
of strong electron donors ((5s=0.38e, blue), and a group of nearly 
neutral atoms (Q_g=-0.02e, green). Due to the undulation, with the 
inequivalent z positions averaging at 0.1254, the Fdd2 structure is 
very similar to a charge density wave (CDW) modulation. The ap- 
pearance of strong charge transfer effects in an element at such high 
pressures is most unexpected, as the resulting Coulomb energy is sub- 
stantial. 

We also found a layered Cmca structure similar to "black phos- 
phorus", which becomes thermodynamic ally stable at about 2 TPa. 
As depicted in the supplementary information, the Cmca structure 
consists of three-fold-coordinated nitrogen atoms and has zig-zag 
layers, with a shortest N-N distance of about 1 . 1 5 A at 2.5 TPa, which 
is much shorter than the N-N separation between the layers of about 
1.56 A. Layered structures also appear in the high-pressure phases of 
other small molecules of first row atoms, e.g., CO |46 1 and CO2 |47 1. 
The nitrogen atoms in the 7?3m structure are six-fold-coordinated 
and form distorted octahedra. Pure four-fold coordinated structures 
often appear in our searches at 3 TPa, but they are less favorable than 
Cmca and Pi/nbm. The crystal structures and charge transfers of 
the newly predicted stable phases are summarized in Table[T| 
Energetics. As can be seen in the enthalpy-pressure relations 
of Fig. [2ja), solid nitrogen undergoes a series of structural phase 
transitions beyond the previously-known framework P2\2\2\ poly- 
meric phase 1251 1261 and recently discovered diamondoid structure 

(743m) /43m ^iZ^^ Cmca ^:ll^f> P4/n&m li^Zf^ 
Fdd2 ^^'^ '^^'^^ RSm ^'^ '^^^^ lAi/amd. This partially ionic 
PA/nbm phase is stable over a wide pressure range, but compres- 
sion to about 7.1 TPa leads to the modulated Fdd2 structure. The 
Fdd2 structure is the most favorable phase from about 7. 1-1 1 .5 TPa, 
whereupon it transforms into a six-fold-coordinated hexagonal i?3m 
phase. As shown in the inset in Fig. [2] nitrogen forms a P4i / amd 
structure similar to Cs-IV |48| at about 30 TPa. 

The volume-pressure relations of the most stable phases are 
shown in Fig.j2|b). The volume decreases by about 2.5% at the tran- 
sition from /43m to Cmca at 2.1 TPa, 0.7% at the transition from 
Cmca to PA/nbm at 2.5 TPa, 0.6% at the transition from PA/nbm 
to Fdd2 at 7.1 TPa, and 0.5% at the Fdd2 to RZm transition at 11.5 
TPa. 



Electronic structures. The electronic band structures and projected 
densities of states (PDOS) of Cmca at 2.0 TPa, P4 /nbm at 3.0 TPa, 
and Fdd2 at 8.0 TPa, are shown in Fig. [5] decomposed into s and p 
components. The P2\2\2\ and /43m structures are semiconducting 
at 2.5 TPa, but the band structures of Cmca, PA/nbm, and Fdd2 
shown in Fig.jsjand R3,m are metallic, see the supplementary mate- 
rial. The insulating state of nitrogen persists to about 2.0 TPa, which 
is considerably higher than in other light elements such as hydrogen, 
carbon and oxygen. f301 

The band structure of Cmca at 2.0 TPa features a partially filled 
Dirac cone between the Y and T points, similar to the one that ap- 
pears in electron doped graphene. |49 1 The projected density of states 
in these metallic structures arises from the s and p orbitals which 
extend over essentially the same energy interval, with both s and p 
orbitals contributing to the electronic density of states at the Fermi 
level and conduction bands, which indicates strong sp hybridization 
in the Cmca, PA/nbm and Fdd2 structures. The densities of states 
at the Fermi levels of PA/nbm and Fdd2 are higher than in Cmca, 
implying that Cmca is less metallic. 

Discussion 

The dynamical stability of the newly predicted stable structures 
was investigated by calculating their phonon dispersion relations. 
As shown in Fig. |4]and in the supplementary material, all of the 
structures are predicted to be mechanically and dynamically stable at 
the specified pressures. The phonon dispersion relation of PA/nbm 
shows steep acoustic branches together with fairly flat optical modes. 
This supports the view that PA/nbm is formed of atomic units: the 
N-N vibrons as well as the breathing mode of the tetrahedra are es- 
sentially dispersionless, and well screened by the electronic cloud. 

On the other hand, the phonon mode of Fdd2 corresponding to 
the amplitude variation of the undulating layers drops to a very low 
frequency of about 190 cm^^. For comparison, the next lowest-lying 
mode (the relative sliding mode) has a much higher frequency of 
about 1300 cm^^, because the shear mode brings atoms of the same 
charge closer together. 

The formation of N2 pairs and N5 tetrahedra in PA/nbm was de- 
duced from the calculated bond lengths and confirmed by the charge 
density plot in Fig. |5|b), where the atoms between the tetrahedron 
center and corners, and between the two atoms within the dimers, 
form strong covalent bonds. The electron localization function (ELF) 
shown in the supplementary material provides additional evidence for 
the formation of tetrahedra and dimers in PA/nbm. 

The formation of charged units from an elementary compound 
has also been observed in 7~Boron |50|. The pressures of interest 
in the present study are much larger than those at which 7~Boron 
has been observed and, besides, P4/7ibm-nitrogen is a very differ- 
ent system. The charge transfer in boron satisfies its tendency to 
form electron deficient icosahedra. Based on the structural richness 
of phosphorus and arsenic one might assume that nitrogen could form 
similarly rich bonding patterns once the triple bond is broken. At 2 
TPa, however, packing efficiency is crucial and one would expect 
structures with high coordination numbers to be formed. For exam- 
ple, the coordination number of Cmca is 3, compared with 6 in the 
/?3m phase. It is therefore reasonable to expect that the average co- 
ordination of nitrogen in the PA/nbm and Fdd2 structures would 
be between 3 to 6. In PA/nbm, if one employs a longer bond cri- 
terion, e.g., 1.30 A, the coordination number for the comers of the 
tetrahedra, their centers and the nitrogen atoms in the chains are 3, 
4, and 6, respectively, as shown in Fig. S5 of the supplementary in- 
formation. In this metallic phase the variation in coordination (or the 
inhomogeneity of the density of ions) leads to charge transfer from 
the highly-coordinated atoms ("N2 pairs") to the lower coordinated 
ones ("tetrahedra" comers). This charge transfer allows the PA/nbm 
stmcture to form a unique metallic all-nitrogen salt. Each atom in the 
ID chains perpendicular to the tetrahedra layers has four long bonds 
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to the comers of the tetrahedra, thus only one electron is left. The ID 
chain with uniform N-N distances is unstable and undergoes a Peierls 
distortion, similar to that in lithium |51 1, and N2 pairs are formed. 

Although Fdd2 also shows clear variation in coordination and 
charge transfer, the atomic arrangement is completely different from 
that in PA/nbm. This phase has undulating layers, formed from 
distorted parallelepipeds, resembling a CDW. CDWs have been ob- 
served previously in chalcogenides under pressure 1521 . In light of 
the charge transfer, the buckling of layers could be viewed not only 
as a symmetry breaking, but also as allowing the positive and negative 
ions to approach one another and reduce the energy. In addition, the 
charge transfer might create atomic sites similar to those in carbon at 
comparable pressures, although the Fdd2 bandstructure shows very 
different behaviour |34|. While carbon expels charge from the 2s 
levels at high compression, eventually forming an electride, the bot- 
tom of the Fdd2 valence band shows a large 2s population. Charge 
depletion from the 2s levels will only arise with the transition to the 
superconducting R3m phase, at almost 12 TPa. Electride structures 
are found for example in carbon 1 34 1 and aluminum |36| at very high 
pressures, but they do not to appear in nitrogen up to at least 100 
TPa. While intuition suggests that close-packed structures should be 
favored under extreme pressures, this is not the case for nitrogen even 
at 100 TPa where, for example, the fee structure is almost 2 eV per 
atom higher in enthalpy than Mi/amd. 

Methods 

Structure search. The ab initio random structure searching (AIRSS) 
method 1431 1441 was used to identify low-enthalpy structures of ni- 
trogen at multi-TPa pressures. Extensive searches were performed at 
selected combinations of 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 
16, 17, 18, 19, 20, 21, 22, 24 nitrogen atoms per cell at 1, 2, 3, 4, 5, 
7, 10, 15, 20, 50 and 100 TPa. More than 47,000 structures in total 
were optimized using DFT. Additional searches were performed us- 
ing initial structures consisting of random packings of N2 dimers and 
N5 tetrahedra, in order to study the partially-ionic structures in more 
detail. No better structure than PA/nbm was found. 

Ab initio calculations. We used the castep plane-wave DFT code 
[531 and an ultrasoft pseudopotential generated with a small core ra- 
dius to avoid overlapping of the core radii at the high-pressures stud- 
ied. The enthalpy-pressure relations of the structures were recalcu- 
lated using very hard Projector Augmented Wave (PAW) pseudopo- 
tentials and the VASP code |54| with a plane- wave basis set cutoff 
energy of 1000 eV. Phonon and electron-phonon coupling calcula- 
tions were performed using DFT perturbation theory |55 | with the 
ABINIT 156 1 and Quantum ESPRESSO codes 157J. Very high 
cutoff energies (1632 eV for ABINIT and 952 eV for Quantum 
ESPRESSO), together with large numbers of k-points, were used to 
obtain convergence of the enthalpy differences between phases to bet- 
ter than 5 x 10^'^ eV per atom. Calculations were performed with 
the different codes to cross-check the transition pressures. 
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Table 1. Structures of the newly predicted stable nitrogen phases. Qm is the IVIulliken charge obtained from ultrasoft pseudopotential 
calculations, and Qs is the Bader charge from PAW calculations. 



Structures 


P (TPa) 


lattice constants 


atomic position 


Qm 


Qb 


Cmca 


2.5 


a=2.090, fe=3.041, c=2.617 


8f (0.0, -0.6034, -0.1843) 












a=90, /3=90, 7=90 








P4:/nbm 


2.5 


a=3.424, 6=3.424, c=2.466 


2d (0.0, 0.5, 0.5) 


0.18 


-0.066 






a=90, 13=90, 7=90 


4g (0.0, 0.0, 0.7376) 


0.11 


0.183 








8m (0.1808, 0.6808, 0.7879) 


-0.10 


-0.075 


Fdd2 


8 


a=5.715, 6=4.590, c=2.404 


16b (-0.8991, 0.9116, -0.8784) 


-0.15 


-0.359 






a=90, /3=90, 7=90 


16b (-0.7687, 0.6347, -0.7130) 


0.14 


0.379 








16b (-0.3059, 0.4036, -0.3075) 


0.01 


-0.020 


R3m 


12 


a=1.613, 6=1.613, c=1.468 


3b (0.0, 0.0, 0.5) 












a=90, /3=90, 7=120 








Mi/amd 


30 


a=0.890, 6=0.890, c=3.452 


4a (0.0, 0.25, 0.875) 









a=90, /3=90, 7=90 




(b) 




Fig. 1. Crystal structures of the newly predicted nitrogen pliases. (a) P4/nbm, 
(b) Fdd2, (c) R3m and (d) /4i /amd. The blue balls in (c) and (d) represent nitrogen 
atoms without obvious charge transfer. The colors in (a) and (b) represent atoms with 
different charges, see the text for more details. The tetrahedra in Pi/nbm, octahedra 
in R3m and trigonal dodecahedron in Mi /amd are depicted in pink, green and cyan, 
respectively. 
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Fig. 2. (a) Enthalpy-pressure relations and (b) volumes. The enthalpies of the 
previously-reported polymeric phases are shown as dashed lines and the solid lines are 
new structures from this work. To make the plot clearer, only the volumes of the ground 
state structures are plotted in (b). The inset in (b) shows the enthalpies of Fdd2 and 
/4i / amd relative to i?3m at very high pressures. 
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Fig. 5. Electron densities of polymeric nitrogen, (a) Cmca at 2 TPa along [100] 
directions, (b) PA/nbm at 2.5 TPa aiong [100] and [110] directions, (c) Fdd2 at 8 TPa 
along [001] directions, (d) R'im at 12 TPa along [110] and [111] directions. 
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